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The density-functional theory (DFT) and density-functional perturbation theory (DFPT) are employed to study 
the pressure dependence of elastic and dynamical properties of zinc-blende ZnS and ZnSe. The calculated elas- 
tic constants and phonon spectra from GPa to 15 GPa are compared with the available experimental data. 
Generally, our calculated values are overestimated with experimental data, but agree well with recent other 
theoretical values. The discrepancies with experimental data are due to the use of local density approximation 
(LDA) and the effect of temperature. In this work, in order to compare with experimental data, we calculated 
and discussed the pressure derivatives of elastic constants, the pressure dependence of dynamical effect charge, 
and mode Gruneisen parameter at r. 
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1. Introduction 

Zinc-blende zinc sulfide (ZnS) and zinc selenium (ZnSe) are wide band II-VI semiconductors. They are 
extremely useful in the manufacture of semiconductor devices. Thorough understanding of the pressure 
dependence of elastic and dynamical properties is the essential prerequisite to a material's synthesis and 
application. The elastic constants of the materials are essential in predicting and understanding the ma- 
terial response, strength, mechanical stability and phase transitions. However, only a few studies have 
been conducted on the mechanical properties of ZnS and ZnSe at elevated pressure as well as pressure 
dependence of their elastic constants. Several research groups have investigated pressure dependence 
of elastic properties of zinc-blende ZnS and ZnSe using different methods. For example, Berlincourt et 
al. (lh using a piezoresonance technique, Lee £2] with ultrasonic echo measurement, and Hodgins Ql us- 
ing Brillouin scattering have determined the elastic constants of ZnS and ZnSe. From measurements at 
different temperatures 80 to 300 K, Lee estimated K values by extrapolation. Theoretically, pseudopo- 
tential plane-wave [PPPW] approach |4], linear muffin-tin orbital [LMTO] method | 5], linear combination- 
atomic-orbitals [LCAO] method |6], full-potential augmented plane wave plus local orbitals (FP-APW+lo) 
method 1 7] etc, have been used to evaluate their mechanical properties under pressure. Several attempts 
have also been made to study the dynamical properties for zinc-blende ZnS and ZnSe. Experimentally, 
the phonon spectra of ZnS and ZnSe (at zero pressure) have been determined respectively at room tem- 
perature using inelastic neutron scattering (INS) by Vagelatos et al. |8] and Hennion et al. |9||. The optical 
phonon modes at the F point have been measured through Raman spectroscopy by Lin et al. 1 10] and 
Hennion et al. @]. Due to an anharmonic temperature shift, the measured frequencies are expected to be 
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somewhat lower than the calculated ones. Several calculations fl ElUlill have been performed to study 
the dynamics of ZnS and ZnSe at zero and high pressures. These calculations can be classified into three 
categories: (1) Model calculations, using the rigid-ion and the adiabatic bond charge models. In these cal- 
culations no softening of phonon modes has been obtained. (2) Ab initio frozen-phonon calculations of 
the frequency of the phonon modes at the T and X points. (3) Density functional perturbation theory 
(DFPT) calculations. Especially Hamdi et al. |4] treated Zn 3d electrons as valence states and using this 
method studied the vibrational and thermal properties of ZnS and ZnSe. Cardona et al. 1 11] using this 
method measured the vibrational properties of ZnS for several isotopic compositions. 

From the above it is clear that there are considerable experimental and theoretical works on ZnS 
and ZnSe. We note that there exist only limited theoretical studies on elastic and dynamical properties 
under pressure. Moveover, The accurate measurement of these quantities is a difficult task due to difficult 
experiment conditions at high pressure. In this paper, we compare the results of experiments with ab 
initio theoretical calculations of the elastic and dynamical properties. We, therefore, think it worthwhile 
to perform these calculations. The rest of the paper is organized as follows: After a brief introduction 
in section 1, the method of calculation and computational details are given in section 2. The simulation 
results for structural, elastic, and dynamical properties are presented and discussed in section 3. Finally, 
the conclusion is given in section 4. 

2. Theoretical method 

The theoretical calculations were performed in the framework of the density functional theory by 
using a local density approximation for the exchange-correlation potential as implemented in the ABINIT 
pac kage jl^Q- Kohn-Sham orbitals are expanded in plane waves with the use of the Troullier-Martins 
11611 pseudopotentials to describe the valence electrons. A 40 hartree cutoff was used and a Monkhorst- 
Pack grid of 12 x 12 x 12 points was used to describe the electronic properties in the Brillouin zone. The 
phonon frequencies were calculated on a 6 x 6 x 6 g-point mesh. These calculating parameters are chosen 
to ensure the total energy error in 0.1 mHa. 

3. Results and discussion 
3.1. Structural properties 

The equilibrium volume of ground state of the zinc-blende phase of ZnS, ZnSe is determined by cal- 
culating the total energy per primitive unit cell as a function of V. The calculated results are shown in 
figure [1] The Murnaghan's equation of state is then used to fit the calculated energy-volume data. The 
obtained structural parameters are compared in table [1] with the available experimental data and other 
theoretical results. The table shows that our calculations underestimate the equilibrium lattice param- 
eter (a) by 1 -5- 2% and overestimate the static bulk modulus (£) by 5 -s- 8% with respect to experimental 
data. Our calculational values are similar to other theoretical values in table [T] The underestimate of 

Table 1. The calculated lattice a, bulk modulus B and their comparison with experiments and available 
theoretical calculations for zinc-blende ZnS and ZnSe. 





ZnS 
a (A) 


B (GPa) 




ZnSe 
a (A) 


B (GPa) 


Present 


5.336 


81.2 


Present 


5.582 


70.8 


FP-APW+lo [7j 


S.342 


89.67 


LAPWF12J 


5.543 


70.0 


PPPW[17] 


5.328 


83.8 


PPPW[4J 


5.543 


70.0 


LMTO [6J 


5.335 


83.7 


LMTO [5J 


5.633 


81.1 


Expt. [8J 


5.412 


75 


Expt. [20] 


5.667 


69.3 
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Figure 1. Calculated total energy per unit cell versus volume for zinc-blende ZnS (left panel) and ZnSe 
(right panel). 



lattice constants compared with experimental values is caused by selection of pseudopotential and LDA 
itself, but these errors are within the acceptable error bars due to the use of LDA 1 19], which reflects the 
reliability of our self-consistent calculations and the pseudopotentials used. On the other hand, the exper- 
imental data are obtained at room temperature and our values neglect thermal expansion. One should 
also notice that the experimental values of bulk modulus are uncertain due to the difficulty of growing a 
high-quality single crystal. 



3.2. Interatomic force constant 

The interatomic force constants (IFC) describing the atomic interactions in a crystalline solid are de- 
fined in real space as lilll 

d 2 E 

Cka.k'P ( fl > W = — — —r . (1) 
01, 01,. a 

ka k'P 

Here, T^ a is the displacement vector of kth atom in the ath primitive cell (with translation vector R a ) along 
a axis. E is the Born-Oppenheimer (BO) total energy surface of the system (electrons plus clamped ions). 
IFC offers a convenient way of storing the information in the dynamical matrix. Furthermore, an appro- 
priate description of the motion of ions in DFPT is necessary. The IFC can be decomposed into an electro- 
static (Ewald) contribution, which is long-ranged, and a "local" contribution which can be attributed to 
covalent bonding. The behavior of the total IFC and of the local contribution as a function of interneigh- 
bor distance is shown in figure[2] Dynamical matrices have been calculated on a 8 x 8 x 8 reciprocal space 
fee grid. Fourier deconvolution on its mesh yields real-space interatomic force constants up to the ninth 
neighbor shell. This procedure is equivalent to calculating real-space force constants using an fee super- 
cell whose linear dimensions are four times larger than the primitive zinc-blende cell, thus containing 
128 atoms. Generally, the decay of local interaction for cation-cation, cation-anion, and anion-anion are 
faster than those of total interaction. At the third neighbor, the local part goes to zero for each species 
pairs of ZnSe, and at the forth neighbor the local part goes to zero for the same species pairs of ZnS, but 
only at the eighth neighbor it goes to zero for any total interaction. The local interatomic force constants 
drop off really rapidly, the first-neighbor force constants being over 5 times as large as any other force 
constants. Since zinc-blende ZnS and ZnSe are polar materials, when an atom is displaced from its origi- 
nal position it creates a dipole. The macroscopic electric field, caused by the long-range character of the 
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Coulomb forces contributes to the longitudinal optical phonons in the long-wavelength (q — 0) limit. The 
real-space interatomic force constants decrease from ZnSe to ZnS, which can be caused by different bond 
lengths and ion strengths. 
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Figure 2. Real-space interatomic force constants for zinc-blende ZnS (upper panel) and ZnSe (lower 
panel). 



3.3. Elastic properties 

The linear elastic constants are formally defined as 



drj rS 



(2) 
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Table 2. Calculated and experimental elastic constants (in GPa) and internal-strain parameter <; for 
zinc-blende ZnS and ZnSe. 





Cll 


Cl2 


C44 


? 


ZnS 


Present 


122 


68 


57 


0.62 


Expt. [1] 


104.0 


65.0 


46.2 




FP-LMTO [5J 


123.7 


62.1 


59.7 


0.651 


FP-APW+lo [7J 


118 


72 


75 


0.715 


PAW[6J 


97.2 


56.4 


64.2 




ZnSe 


Present 


91.2 


58.2 


42 


0.61 


Expt. [2] 


85.9 


50.6 


40.6 




Expt. [22J 


82.8 


46.2 


41.2 




FP-LMTO [5J 


95.9 


53.6 


48.9 


0.63 


FP-APW+lo [7J 


94 


61 


64 


0.746 


PPPW[4J 


91.3 


56.3 


38.3 


0.734 



where a and r\ denote the externally applied stress and strain tensors, respectively. In the case of cubic 
crystals, there are only three independent elastic constants, C\\ = c X x,xx> c i2 - c X x,yy, C44 = Cy Z ,y Z - The 
en, c\2 and unrelaxed (i.e., in the absence of any internal displacement) C44 elastic constants have been 
calculated by computing the second derivative of BO energy surface with respect to the strain, namely 

1 d 2 E tot (V P ) 

Cap CP) = 77 T—^ . (3) 

r Vp OE a OEp 

where c a p [P) is the pressure dependence of the elastic constants, E tot ( Vp) is the total energy per unit cell, 
Vp is the unit cell volume at a given pressure P. The calculations of pressure dependence of the elastic 
constants are performed in two steps. Firstly, we calculate the total energy of the bulk crystal as a func- 
tion of the unit cell volume. Then, using the definition of pressure, P = -dE tot /dV, we calculate the unit 
cell volume corresponding to a certain value of the external pressure P. In the second step, the unit cell 
at certain pressure is subjected to test distortions. To determine the pressure-dependent elastic constant, 
the deformation energy has been computed for a series of test displacements in the range of -1% to 1% 
and have been fitted by the second-order polynomials to the expressions from the elasticity theory. The 
theoretical values of elastic constants and of the internal-strain parameter of zinc-blende ZnS and ZnSe 
are summarized in table [2] The calculated elastic constants are overestimated by about 5% to 15% with 
experimental data obtained from Brillouin scattering measurements. These errors with experimental 
data are within the acceptable error bars due to the use of LDA. The experimental data are obtained at 
room temperature mostly, which is easy to understand from the fact that increasing temperature leads 
to an increase of lattice constants, and finally leads to a decrease of elastic constants. Good agreement 
is observed with other recent theoretical values. To the best of our knowledge, no experimental data for 
the internal-strain parameter of both materials are available. Although our calculated internal-strain pa- 
rameters of zinc-blende ZnS and ZnSe are in good agreement with recent FP-LMTO |5] and FP-APW+lo tJ 
calculations, their absolute values are slightly smaller than the previously calculated ones. 

In figure [3] we present the pressure dependence of elastic constants, c\\{P), C12CP) and C44OP), ob- 
tained for zinc-blende ZnS and ZnSe in the range of hydrostatic pressure -=- 15 GPa. The figure displays 
the calculated points (circles), and solid lines are quadratic fits to the calculated results. One can notice 
that the elastic constants c\\ and C12 increase with pressure more significantly than C44. Good agreement 
is found for the pressure derivatives of cu with measurement up to 0.6 GPa for ZnSe | 2]. At P - 0, we find 
for ZnSe: dcn/dP = 4.33, close to the experimental value 4.44, dci2/d_P = 4.62 also in good agreement 
with the measured value 4.93. For dc44/d_P we found 0.79 which is somewhat larger than the measured 
0.43. The discrepancy is due to temperature effect. Indeed, the experimental data of reference |2j are 
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Figure 3. Pressure dependence of elastic constants Cfj for zinc-blende ZnS (left panel) and ZnSe (right 
panel). Open diamonds, asterisks and squares are taken from experimental data from reference (J, 
and 1 22], respectively. The solid lines are quadratic fits to the calculated results. 



obtained at 295 K. It is found theoretically by ab initio calculation that the elastic constants decrease with 
temperature increase |4]. The calculated hydrostatic pressure coefficients at P = for ZnS are 3.88, 4.76 
and 1 for C\i, Cyi and C44, respectively. To our knowledge, there is no experimental data available for 
dcij/dP of ZnS. 



3.4. Dynamical properties 

The calculated phonon dispersion curves at P = GPa and P = 9 GPa for zinc-blende ZnS and ZnSe 
are displayed in figure |U The vibrational frequencies were determined at several volumes within the 
linear response framework. There is no gap between the acoustical and optical phonon branches for 
ZnSe. The overlap is caused by the nearly identical masses of Zn and Se atoms. In figure [4] our result 
is compared with experimental data at ambient pressure from reference 

SEE! 

. In particular, for 

ZnSe, the transverse optical (TO) and longitudinal (LO) phonon modes at P = GPa at zone center are 
found to be 216 cm -1 and 253 cm -1 , and those phonon modes have been reported with frequencies of 
213 cm -1 and 253 cm -1 in inelastic neutron scattering |9], but our TO at Y is larger than the value in 
Raman scattering 1 10] . For ZnS, the calculated phonon dispersion at P = GPa is overestimated with 
experimental data 1 23] except for the transverse acoustical (TA) phonon modes at L in Raman scattering. 
These differences may be enhanced by the fact that the LDA value of ag is smaller than the experimental 
one. Our calculated phonon frequencies are in good agreement with other theoretical results for both 
ZnS (lllQ and ZnSe (HIH. 

The pressure dependence of LO, TO and LO-TO splitting are shown in figure [5] In a polar lattice, the 
splitting of the optical phonon modes is determined by two parameters, i.e., Born's dynamical effective 
charge of the lattice ions and the screening of the Coulomb interaction, which depends on the electronic 
part of the dielectric constant in the phonon frequency regime. In this work, the phonon frequencies shift 
to higher energies and LO-TO splitting decreases with pressure. The change of the dynamical effective 
charge under pressure can be determined from the frequencies of the optical phonon using equation 1 26] 



Z* 2 = £n£oo^M( w LO ~ W TC>) ' 



(4) 



where £n is the vacuum permittivity, \i is the reduced mass of an anion-cation pair. V is the available 
volume per pair, and to is angular mode frequency. The calculate dynamical effective charge as a function 
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Figure 4. Calculated phonon dispersions for zinc-blende ZnS (upper panel) and ZnSe (lower panel) at 
zero (solid lines) and 9 (dotted lines) GPa pressures. Experimental data are denoted by open circles and 
asterisks from reference (2^] for ZnS, reference 0] (open circles), (asterisks) for ZnSe. 

of pressure are shown in figure[6] The calculated effective charge decreases slightly with rising pressure. 
This finding indicates that a decrease of the Born's dynamical effective charge at high pressure is the sign 
of strong covalent bonding and the related overall increase of direct optical gaps with pressure. 

The pressure dependence of phonon frequency is usually expressed in terms of the mode-Gnineisen 



Table 3. Mode-Griineisen parameters at r point for ZnS and ZnSe. 





ZnS 
LO 


TO 




ZnSe 
LO 


TO 


Present 


1.12 


1.38 


Present 


1.04 


1.57 


Expt. 


0.90 a 


1.27° 


Expt. 


0.90 fc , 0.85 c 


\A h , 1.52 c 


Theo. 


1.12 & 


1.66 & 


Theo. 


1.18 d , 1.19 & 


1.47 d , 1.66 & 



reference [27], b reference [4], c reference (28|], d reference [29]. 
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parameter 

Yj.q = -{d\nwj, q /d\nV) = {B 1 w j , q ) {dh> j , q / dP) . (5) 

Available theoretical and experimental data of mode-Griineisen parameters of the optical T mode 
are collected in table [3] Our ab initio values of optical T mode-Griineisen parameters seem to overesti- 
mate the experimental value. Assuming anharmonic effects other than those of thermal expansion to be 
negligeable, we define an effective mode-Griineisen parameters as 

1 din a din a 

f- , aj — . (6) 

3ar dT dT 

From temperature dependence data and using the coefficient of linear thermal expansion, aj = 6.8 x 
10 -6 K -i fQr ZnS Q and ar _ 6 34 x 10 -6 K -i for ZnSe at 3QQ K> we est i mate f TQ = L31 f or ZnS and 

7to = 1-53 for ZnSe at room temperature, slightly smaller than our yxo- 



400 





30 1 1 1 1 

6 10 16 

Pressure (GPa) 

Figure 5. Pressure dependence of the zone-center optical phonon frequencies (upper panel) and differ- 
ence of the optical phonon frequencies (lower panel) for zinc-blende ZnS (open circles) and ZnSe (aster- 
isks). The solid lines are quadratic fits to the calculated results. 
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Figure 6. The pressure dependence of dynamical effective charge for zinc-blende ZnS and ZnSe. The solid 
lines are quadratic fits to the calculated results. 

4. Conclusion 

We present the results of pressure dependence of elastic and dynamical properties of zinc-blende 
ZnS and ZnSe. The considered P ranges are + 15 GPa. The calculated elastic constants c\\ and c\2 in- 
crease with pressure more significantly than C44. The calculated values are overestimated by about 5% to 
15% with experimental data obtained from Brillouin scattering measurements, but there is good agree- 
ment with other theoretical values. The errors with experimental data are due to the use of LDA, and 
the increasing temperature in the experiment leads to a decrease in the elastic constants because the ex- 
perimental data are mostly obtained at room temperature. In this work, the calculated TO at f is larger 
than the value in Raman scattering for ZnS, and the calculated phonon dispersion at P — GPa is over- 
estimated with experimental data except for the transverse acoustical (TA) phonon modes at L in Raman 
scattering. By defining effective mode-Griineisen parameters 7 and using the coefficient of linear thermal 
expansion, we see that fro for ZnS and ZnSe at room temperature is slightly smaller than our jto- 
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3a/ie>i<HicTb Bifl Ti/iCKy npy>i<Hi/ix i fli/maMmm/ix B/iacTi/iBOCTefi 
141/1 hkoboi o6i\/iaHKi/i ZnS i ZnSe 3 nepujonpi/mminHHX 
po3paxyHKiB 

X.M. Bam™, A>k.M. KecR CM. XyamE L1..M. XyanP 

(PaKy^bTeT cf>i3HKH Ta e/ieicrpoHiKn, yHiBepci/ueT KciaHaH, HeHb4>K0y 423000, KHP 
tt>i3M4Hnti 4>iKy^bTeT, LleHTpa/ibHO-niBfleiHHufi yHiBepci/ueT, HaHiiia 41 0083, KHP 

Teopia cfiyHKuiOHa/iy rycTWHW (DFT) i Teopm 36ypeHb cfiyHKuioi-iayiy rycTWHW (DFPT) 3acrocoByeTbCfl flo BUBHe- 
hhh 3a/ie>KHOCTi Bifl TUCKy npyxHwx i fli/maMiMHUx B/iacTHBOCTew uuhkoboi o6MaHKH ZnS i ZnSe. Po3paxoBaHi 
npy>KHi cra/ii i cfiOHOHHi cneKTpn Bifl GPa flo 15 GPa nopiBHioiOTbCfl 3 HanBHUMW eKcnepuMeHTa/ibHUMW fla- 
humh. B3ara^bHOMy, Haiui po3paxoBaHi 3Ha4eHHfl e nepeoi^iHeHi b nopiBHflHHi 3 eKcnepuMeHTayibHUMn, a/ie 
y3roflxyK)TbCfl flo6pe 3 imui/iMW TeopeTH4HHMn 3Ha4eHHflMn. Po36i>KHOCTi 3 eKcnepuMeHTayibHUMW flaHWMn e 
nepe3 BMKopucTaHHfl Ha6/inxeHHfl /lOKa/ibHOi rycTHHH (LDA) i Bn/11/iBy TeMnepaTypn. fl,nn nopiBHflHHfl 3 eKcne- 
puMeHTa/ibHUMn flaHUMW, Mi/i o64waiK)eMO i o6roBopioeMO noxiflHi 3a tuckom Bifl npyxm/ix da/inx, 3a/ie>KHiCTb 
Bifl TUCKy flUHaMiHHOro ecfieKTHBHoro 3apafly i napaMeip rptOHaifaeHa npw T. 

lOnoHOBi c/ioBa: ab initio, crpyKTypa, MixamMHa cwnoBa cra/ia, tuck, npyxHicrb, umaMiKa 
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